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Abstract. Upcoming wide-area surveys in the submillimetre regime will allow the construc- 
tion of complete galaxy cluster samples through their thermal Sunyaev-Zel'dovich effect. We 
propose an analytic method to predict the number of gravitationally lensed giant arcs produced 
by the cluster sample expected to be detectable for the Planck satellite. The statistics of lensed 
arcs can be used to constrain cosmological parameters. This computation implies the choice of 
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Q-r an analytic cluster model which has to reproduce the lensing properties of realistic galaxy clus- 
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ters. Comparing the results of analytic computations and numerical ray-tracing simulations, we 
show that spherical models with the Navarro-Frenk- White density profile fail by far to produce 
as many giant arcs as realistic galaxy clusters. This confirms the great importance of asymme- 
tries and substructures in the lensing matter distribution. By elliptically distorting the lensing 
potential of spherical NFW density profiles, we show that the discrepancy between analytic and 
J> ' numerical computations can almost completely be removed. 
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1 Introduction 

Several studies demonstrated that arc statistics is a very sensitive tool to constrain cosmological mod- 
els. For example, using ray-tracing simulations and numerically simulated galaxy clusters, Bartelmann 
et al. jjj showed that the expected number of giant arcs (defined as arcs with length-to-width ratio 
larger than 10 and B magnitude less than 21.5) changes by orders of magnitude from low density to 
high density universes, depending also on the value of the cosmological constant. This result, with 
the exception of the dependency on the cosmological constant, was confirmed by several other authors 
Hi H]> w ho used spherical analytic and somewhat idealistic models to describe the cluster lenses and 
predict the number of giant arcs as a function of the cosmological parameters. Therefore, the com- 
parison of the expectation from such theoretical studies to the observations of strong lensing events 
in galaxy clusters can be used to constrain the cosmological model. 

In order to make such a comparison, a complete sample of galaxy clusters could be selected through 
the thermal Sunyaev-Zel'dovich effect (hereafter SZ effect). This will be possible in the near future 
thanks to the upcoming full-sky microwave surveys, like that planned for the Planck satellite mission. 
Although other selection methods could be used for this purpose (i.e. through the X-ray emission from 
the intra-cluster gas), the SZ selection criterion is much less sensitive to biases due to the dynamical 
processes in the cluster centre, like cooling flows or merging events. 

2 Lensing Sunyaev-Zel'dovich clusters 

Upcoming full-sky microwave surveys will detect of order of 10 4 galaxy clusters through their thermal 
SZ effect. This effect is determined by the Compton-y parameter 

V{6) = I dl T(9, 1) n e {6, 1) , (1) 

where T is the temperature of the ionized intra-cluster gas and n e is the electron number density. The 
clusters detectable through their thermal SZ effect are those for which 



Y = J d%(9) > Y m 



where the integral is performed over the source area. The lower limit Y m { n depends on the angular 
resolution and the temperature sensitivity of the detector. For example, for the Planck satellite, Y m { n 
is expected to be ~ 10 -4 arcmin 2 . 

As explained in more detail by Bartelmann ||, many SZ clusters will also be efficient weak gravi- 
tational lenses. Their weak lensing effects can be quantified using the aperture mass [O], 

mm = J d 2 Mm(m) , (3) 

which is the integral within a circular aperture with radius 9 over the convergence k, weighted by a 
suitable filter function satisfying the condition 

/ d 2 MU(-d) = . (4) 
Jo 

The dispersion of M ap due to the finite number density n g of randomly distributed background 
galaxies and their intrinsic ellipticity dispersion a € was computed by Schneider [fTH| , 

,„(,) = 0.016 ( a V' /2 («) (^J—Y 1 . (5) 
y ' V 30 arcmin 2 / \0.2j \1 arcmin J y ' 

Clusters produce a significant weak-lensing effect if the signal-to- noise ratio S = M ap ((9)/<7M(#) is 
larger than a minimal value S m \ n . Assuming S m \a = 5, Bartelmann Q estimated that more than 70% 
of the original SZ cluster sample expected from Planck will be efficient at least for weak lensing, and 
many of them will also show strong lensing features. 



3 Number of arcs 



The number of arcs which is expected to be observed in the selected cluster sample can be estimated 
using the following information: 

• both conditions, that clusters must be detectable through the thermal SZ effect and be efficient 
weak gravitational lenses, give a minimal cluster mass M m i n (zL) at any given redshift zl §; 

• the number density of dark matter haloes npg(M, z\) can be modeled using the Press & Schechter 
[10 1 mass function or its generalizations, like the Sheth & Tormen [12| mass function; 



the lensing cross section for arcs with a minimal length-to- width ratio L/W (approximately 
corresponding to a minimal magnification /j>) a(L/W, \x\ M, zl, £s) can be computed for any 
galaxy cluster with mass M at redshift and for any source redshift zs once the analytic model 
describing the cluster lens is chosen; 



the source galaxy redshift distribution can be modeled as 

p(z s ) oc z$exp(-zg) . 

where ^ = 1.5 @; 



(6) 



finally, the differential number density of sources with a given observed flux F can be assumed 
to be a power law: 

dns „ _ (7) 
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where a ~ 1.8 in the R-band (IT 
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Using these ingredients, the optical depth for the formation of arcs with length-to- width ratio 
larger than L/W is given by 



t(L/W,h;zs) 
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Figure 1: Comparison between numerical and spherical analytic models: the dashed line indicates 
the cross section for arcs with length-to-width ratio > 7.5 of NFW spheres as a function of the lens 
mass; filled circles indicate the cross sections of the five numerical clusters used in the ray-tracing 
simulations. 



where D$ is angular diameter distance to the sources at redshift z$, and dV/dz is the proper volume 
per unit redshift. 

Finally, the number of arcs with L/W > (L/W) m i n and magnitude R < i? max (corresponding to 
the minimal observed flux Fq) is given by 

N = T T r(L/W^;z s ) ( f° ^fr— ) p(z S )d^z S • (9) 



4 Which analytic model? 

The crucial point in this computation is the choice of an analytic model for describing the lenses. 
Several authors [||, ||] modeled them as isothermal spheres in order to measure their efficiency to 
produce giant arcs. In this study, we use the much more realistic Navarro-Frenk- White (hereafter 
NFW) density profile, and we compare the strong lensing properties of this analytic model with those 
of realistic galaxy clusters obtained from TV-body simulations. The NFW profile has strong-lensing 
properties which differ substantially from the singular isothermal profile || [|. 

For this purpose, we use a sample of five clusters obtained from dark matter iV-body simulations, 
kindly made available by the GIF collaboration. These are placed at redshift z^ = 0.27 and have virial 
masses ranging between ~ 3 x 10 14 M e /h and ~ 10 15 M e /h. The cosmological model is the ACDM 
model (Om = 0.3, = 0.7, h = 0.7, erg = 0.8). A complete description of the numerical cluster 
sample can be found in Bartelmann et al. |l|] . We use that sample to perform ray-tracing simulations 
[0j §]> placing a large number of elliptical sources at redshift z$ = 1. Finally, we measure the lensing 
cross sections for arcs with length-to-width ratio larger than 7.5. 

The results are shown in Fig. |l[ The plot displays the numerically measured cross sections (filled 
circles), obtained by averaging the values corresponding to three independent projections of each 
cluster, as a function of the virial mass of the lens. The error bars show their uncertainties. The dashed 
line indicates the analytically computed lensing cross section of NFW spheres placed at the same 
redshift of the simulated clusters. The numerical cross sections are systematically and substantially 
larger than the analytic ones by roughly two orders of magnitude. This clearly shows that spherically 
symmetric models, even with a very realistic density profile, fail to reproduce the strong-lensing 
properties of numerically modeled galaxy clusters. 

A possible explanation for this discrepancy is the lack of asymmetries and substructures in the 
matter distribution of the spherical models. In fact, the importance of both these factors for the strong 
lensing efficiency of galaxy clusters was demonstrated by Bartelmann et al. [g]. 



Figure 2: Deflection-angle maps obtained assuming a spherically symmetric NFW lensing potential 
(left panel), and after adding an ellipticity of e = 0.3 (right panel), as explained in the text. 



5 Elliptical NFW lensing potential 

The lack of asymmetries and substructures in the matter distribution could be compensated by adding 
ellipticity to the lensing potential of the NFW spheres. 
The spherical NFW lensing potential is given by 



tf(a?) oc ( - In 2 | - 2arctanh 2 y'i-^j , (10) 

where x = r /r s is the distance from the lens centre r divided by the scale radius r s of the dark- matter 
halo. 

We introduce the ellipticity e = 1 — 6/aby substituting 

(11) 




The components of the deflection angle are then 

dtf x x d^> x 2 (l-e) 

ox\ (1 — e)p 0x2 p 

where a(p) is the unperturbed (i.e. spherical) deflection angle at the distance p from the lens centre. 

Using ([l2]), we produce deflection angle maps like those displayed in Fig. |2l The left panel of the 
figure shows the unperturbed deflection angle map for a lens of 7.5 x 10 14 M Q //i at redshift = 0.3. 
The right panel shows the map obtained by adding an ellipticity of e = 0.3. Using these deflection-angle 
maps in the ray-tracing simulations, we measure the lensing cross sections of the elliptical models. 
These are displayed in Fig. || as a function of the ellipticity added. The results show that an ellipticity 
of e = 0.3-0.4 could suffice to remove the discrepancy between the numerical and the analytic models. 
Therefore, the usage of elliptical models is more appropriate for computing the number of arcs with a 
minimal length-to-width ratio produced by a sample of galaxy clusters via gravitational lensing. 



6 Conclusions 



Upcoming full-sky microwave surveys will detect several thousands of galaxy clusters through their 
thermal SZ effect. Many of these clusters will also be efficient gravitational lenses. 
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Figure 3: Lensing cross section for arcs with length-to-width ratio larger than 7.5 as a function of the 
ellipticity added to the lensing potential. The lens mass is 7 x 10 M@/h, and its redshift is zl = 0.3. 




We can predict how many giant arcs will be produced by this cluster sample by using analytic 
models to describe the lenses. We verified that spherical models, even with realistic profiles, do not 
reproduce the lensing properties of numerically modeled galaxy clusters: We need to include at least 
ellipticity in order to give a correct description of the cluster lenses. 

We now plan to calibrate the analytic lensing cross sections for the asymmetric NFW halo profile 
with numerical cluster models. This can be done by measuring the ellipticity of the lensing potential 
of the simulated clusters by fitting their deflection angle maps to those of perturbed analytic models. 
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